Two-parameters electromagnetic hysteresis of a subwavelength nonlinear e-near-zero 

slab 
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We consider propagation of transverse magnetic electromagnetic waves through a nonlinear meta- 
material slab, of subwavelength thickness and very small average dielectric permittivity, made up 
of alternating metal and dielectric nonlinear layers. Exploiting the effective-medium approach, we 
prove that the output field intensity is a multi-valued function of both the input field intensity and 
incidence angle so that the transmissivity exhibits a novel kind of multi-stability and a pronounced 
two-parameter hysteresis behavior. The predicted two-parameter hysteresis physically stems from 
the fact that, since the linear and nonlinear contributions to the overall dielectric response can 
be comparable, the slab can host more than one electromagnetic configuration compatible with 
the incident field and the boundary conditions, regardless of the slab thickness. We investigate 
slab transmission also through full wave simulations proving the robustness of the two-parameter 
hysteresis beyond the effective-medium approximation. 

PACS numbers: 78.67.Pt, 42.65.Tg 



Hysteresis is probably one of the most intriguing fea- 
ture a nonlinear system can exhibit both theoretically 
and for its potential applications for the designing of 
memory devices. Optical bistability and related hystere- 
sis have attracted a large research interest in the last 
decades [ll Q since they are the basic ingredients for de- 
vising optical memories, logic gates and optical comput- 
ing devices Q. The advent of metamaterials, with the 
related possibility of engineering the medium electromag- 
netic response [4j , has triggered a renewal of interest to- 
wards optical bistability Q. One-dimensional photonic 
crystals consisting of alternating layers of positive-index 
and negative-index materials index [7] allow, in the pres- 
ence of a nonlinearity, bistable switching and tunable 
nonlinear transmission omnidirectional bistability H 
and bistability at very low values of input intensity [9|. 
Optical multistability have also been predicted in the 
presence of metal-dielectric multilayer structures [lo| and 
it has been shown that, since for these structures the sign 
of the effective dielectric constant depends on the field in- 
tensity, steplike transmission of light can be achieved [ll| . 
On the other hand, media with very small permittivity 
have been shown to support the regime where linear and 
nonlinear polarizations are comparable [l^ and this has 
led to the prediction of a pronounced transmissivity di- 
rectional hysteresis exhibited by nonlinear metamaterial 
slabs with very small linear permittivity [l3{ . 

In this Letter, we investigate the electromagnetic 
transmission of transverse magnetic (TM) waves through 
a metamaterial slab of subwavelength thickness made up 
of alternating metal and dielectric nonlinear Kerr layers 
stacked along the direction orthogonal to the thickness. 
The layers dielectric properties and filling fractions are 
chosen so as to yield a negative and very small effective 
linear permittivity. Exploiting the effective medium ap- 




FIG. 1: (Color online) (a) Geometry of the metamaterial slab, 
(b) Geometry of the TM incident (i), reflected (r) and trans- 
mitted (t) plane waves together with the TM field within the 
slab. 



proach, we numerically solve Maxwell equations for the 
problem of reflection and transmission of an inclined inci- 
dent TM plane wave and we show that the slab transmis- 
sivity is a multi-valued function of both the input field 
intensity and incidence angle. The predicted multistabil- 
ity physically stems from the fact that, since linear and 



nonlinear dielectric polarizations can be comparable [12 1 
(extreme nonlinear regime), the overall nonlinear dielec- 
tric response supports more than one field vector con- 
figuration compatible with the incoming plane wave and 
the boundary conditions, regardless of the slab thickness. 
Therefore, the proposed multistability fundamentally dif- 
fers from standard optical bistability where the multiple 
field configurations are produced by a nonlinear feedback 
process. In order to support our results beyond the effec- 
tive medium approximation, we perform full wave simu- 
lations whose results are in good agreement with those of 
the effective medium approach thus proving the robust- 
ness of the predicted phenomenology. 

Consider a TM monochromatic field (with time de- 
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pcndence exp(— iwt)) E = Ex{x,z)ex + Ez{x,z)ez and 
H = Hy{x,z)ey propagating through the metamate- 
rial slab, embedded in vacuum, reported in Fig. 1(a) of 
sub- wavelength thickness (along the z-axis) L = 2/kQ 
(ko ~ w/c) and consisting of alternating, along the y-axis, 
two nonlinear isotropic and non-magnetic Kerr layers of 
thicknesses di and c?2- We here consider the situation 
where the layers alternation period is much smaller than 
the vacuum wavelength (i.e. di -f ^2 ^ 27r/fco) so that, 
exploiting the effective-medium approach, the slab be- 
haves as an homogeneous medium characterized by the 
dielectric response (constitutive relations) 



Dx = eo je-B^; + X 



(E.E*)i?, + i(E.E)i?: 

(E.E*)£;, + i(E.E) e; 



(1) 



where the effective dielectric permittivity e — fei -I- (1 — 
/)e2 and effective nonlinear susceptibility x = fxi + (1 ~ 
/)X2 are the averages of their layers counterparts (here 
f_= di/{di + d2) is the volume filling fraction of layer 1) 
|12l |. Therefore, employing both nonlinear negative and 
positive dielectrics (i?e(ei) < and Re{e2) > 0), condi- 
tions can be found so that the effective dielectric permit- 
tivity e has a negative real part (with a very small imag- 
inary part) and very close to zero, and x > (see below 
for an explicit example) which is the situation we focus on 
in this Letter. A TM plane wave (E^, H^) incoming from 
vacuum (z < 0) and impinging on the slab at z = with 
incidence angle 9, produces a TM reflected plane wave 
(Er,Hr) for z < and a TM transmitted plane wave 
(Ei, Ht) for z > L, as reported in Fig. 1(b). The fields of 
the three waves are Eg = Es{ey x kg/fco) exp (ik^ -r), 
Hs = Es a/ eo/noBy exp (ik^ • r) where s = {i,r,t), 
Ei, Er and Et are the amplitudes of the three 
waves and k^ = kt = fco (sin Oe^ + cos Oe^ ), k,. = 
kQ{sin6ex —cosBez) are the wave vectors of the three 
waves. The TM electromagnetic field within the slab 
is of the form E ~ [Ex{z)ex + Ez{z)ez]exp(ikoxsm9), 
H = [Hy{z)ey] exp(ikQXshi6) so that. Maxwell equations 
V X E = iLdfioH, V X H = —iujT) yields the system 



dEx 
dz 

dHy 

dz 

(fco sin 0)Hy 



z(fco sm6)Ez + iiOfioHy, 
iusDx, 



(2) 



where and are given by Eqs.(IT]). At the slab in- 
put z — and output z = L faces we require the con- 
tinuity of the tangential components of the electric and 
magnetic fields. In order to evaluate the slab transmis- 
sivity T = we specify the amplitude of the 

transmitted field Et at the output face z = L and suc- 
cessively solve Eqs.Q all the way to the input face at 
z = 0. For our numerical evaluations we have chosen a 
slab with e = —0.05 and x > and we have discarded 
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FIG. 2: (Color online) (a) Nonlinear (solid line) and linear 
(dashed line) transmissivity T as a function of the normalized 
input field amplitude |xl^'''^l^i| '''t the fixed incident angle 
61 = 0.17 rad for a slab with e = -0.05 and x > 0. (b) Surface 
Ixl^'^^^z obtained by Eq.®, for a slab with e = —0.05 and 

x>o. 
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FIG. 3: (Color online) (a) Nonlinear (solid line) and linear 
(dashed line) transmissivity T as a function of the incident 
angle 9 for the normalized input field amplitude |xr''^l^i| = 
0.05 for a slab with e = -0.05 and x > 0. (b) Surface \x\^^'^A^ 
obtained by Eq.(|3}, for a slab with e = —0.05 and X > 0. 



the solutions of Eqs.([2]) for which the profile of E^ is dis- 
continuous within the slab < z < L (a possibility stem- 
ming from the fact that the third of Eqs.Q is a nonlin- 
ear algebraic equation and not a differential equation for 
Ez). In Fig. 2(a) we report the evaluated transmissivity 
T (solid line) as a function of the normalized input field 
amplitude jxl^^^l-E",! at the fixed incidence angle 6 = 0.17 
rad. In order to evaluate the transmissivity for a given 
input field amplitude as a function of incidence angle 9 
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we have first evaluated a number of amplitudes Ei for 
various Et and 6 and subsequently collected all the con- 
figurations corresponding to the same \Ei\. In Fig. 3(a) 
we report the resulting slab transmissivity T as a func- 
tion of 6* for \x\^^^\E,\ = 0.05 (solid line). In Fig.2(a) and 
3(a) we have also reported, for comparison purposes, the 
linear (x = 0) slab transmissivity (dashed line) obtained 
in the same electromagnetic excitation conditions. Note 
that both the transmissivities (solid lines) of Fig. 2(a) and 
3(a) are multi-valued functions or, more precisely, there 
are ranges of both the normalized input field amplitude 
(0.025 < \x\'^^^\E,\ < 0.035) and the incidence angle 
(—0.1 rad < 9 < 0.1 rad) at which the transmissivity 
is threefold. 

In order to discuss the origin of the predicted multi- 
stability it is worth noting that the first two of Eqs.® 
are differential equations describing the spatial dynam- 
ics of the fields E^ and Hy whereas the third of Eqs.([2]) 
is a nonlinear algebraic (cubic) equation yielding a field 
constraint which allows to evaluate E^ (at least locally) 
for Ex and Hy prescribed. Therefore, it is evident 
that, in general, there are three values of Ez compati- 
ble with given E^ and Hy so that, from a physical point 
of view, the overall nonlinear response of Eqs.([T]) can 
generally support different electromagnetic field config- 
urations for a given input intensity. In order to clarify 
this point and to discuss the impact of the boundary 
conditions, note that, from the third of Eqs.(l2|), it is ev- 
ident that the continuity of Hy across the slab bound- 
aries is equivalent to the continuity of D^, so that, af- 
ter defining the real amplitude through the relation 
Ez {L) = Az {Et l\Et\) exp(2i cos 0) and exploiting the sec- 
ond of Eqs.lH)) together with the continuity of E^, the 
matching condition Dz{L^) = Dz{L^) yields 



tl^cos^i 



Az + \Et\sine = (3) 



which is a relation that, for a given Et and 9 yields the 
possible values of Az- Note that, in the linear regime 
(X = 0) Eq.© yields Az = -\Et\ s\n9/e i.e. the bound- 
ary conditions allow only one possible value of Az. The 
possible values of IxT^^^z obtained by Eq.Q belong to 
a surface which is plotted in Figs. 2(b) and 3(b), in the 
case e = —0.05 and x > 0, and we note that such sur- 
face is folded so that there are regions where, for a given 
pair Ixl^^^l-Etl and 9, there are more than one acces- 
sible value of \x\^^'^Az. Therefore the multi-valuedness 
of the slab transmissivity corresponds to different fields 
with the same input intensity and incidence angle hav- 
ing different output longitudinal component Az belong- 
ing to different sheets of the folded surface. Note that. 



for 



0, Eq.Q admits the solutions A 



(0) 



and 



}(\Et?) 



so that the 



lxr/24*^ = ±V-sign(x) (2e/3- 
surface can be folded only if sign(ex) = — 1 and its upper 
and lower sheets are such that x^z ~ kl- This implies 
that, as opposite to the situation we are considering in 
this Letter where |e| <IC 1, in standard materials where e 
is of the order of unity, the required nonlinearity would 




FIG. 4: (Color online) Transmissivity hysteresis (solid lines) 
obtained by varying 6 and X^^'^\Ei\ along a closed path (re- 
ported as a dashed line). 



be so large to effectively forbid the folding of the surface 
(and the consequent feasibility of the directional multi- 
stability). 

The predicted multistability is accompanied by a two- 
parameters transmissivity hysteresis. Note that each of 
the points of the curves (solid lines) of Figs. 2(a) and 3(a) 
corresponds to an electromagnetic configuration within 
the slab and hence these curves have an image on the 
surface IxP^'^^z which are reported in Figs. 2(b) and 
3(c) as solid lines as well. Suppose to switch on the in- 
cident plane wave with 9 = 0.17 rad and with an input 
intensity very small so that the state A of Fig. 2(a) and 
2(b) is excited since before the illumination the slab was 
not polarized. By increasing the input intensity and hold- 
ing the incidence angle fixed, the transmissivity follows 
the curve joining the point A and B of Fig. 2(a) while the 
electromagnetic state continuously varies as in Fig. 2(b). 
From the state B, if the intensity is further increased, 
it is evident that the electromagnetic state undergoes a 
sudden jump to the state C (see Fig. 2(b)) on the lower 
surface sheet since there is no allowed state continuously 
joined to B. As a consequence, the transmissivity under- 
goes a sudden jump to a higher value (see Fig. 2(a)). If 
now, starting from the state C, the input intensity is de- 
creased, the electromagnetic state varies along the curve 
from C to D of Fig. 2(b) whereas the transmissivity in 
Fig. 2 (a) assume the values along the curve from C to D 
which are different from those attained along the forward 
path. At the state D, if the intensity is further decreased, 
the state undergoes a jump to the state E of Fig. 2(b), be- 
longing to the central surface sheet, since the states from 
D to F of Fig. 2(b) have, according to Fig. 2(a), higher in- 
put intensities. An analogous but different hysteretic be- 
havior is obtained if the input intensity is held fixed and 
the incidence angle is varied, as depicted in Fig. 3. Since 
the presented multistability depends on two parameters, 
one can vary both the incidence angle and intensity on 
a prescribed path to obtain nontrivial hysteresis behav- 
iors, as reported in Fig. 4 where 9 and x^^^l^-il are varied 
along a closed path. 

In order to discuss the feasibility of the above pre- 
dicted two-parameters multistability we have performed 
3D full-wave finite-element simulations [l3| for evaluat- 
ing the transmissivity of the slab reported in Fig. 1(a) in 
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FIG. 5: (Color online) Comparison between the slab transmis- 
sivities evaluated through full-wave simulations (dotted lines) 
and those of Fig. 2(a) and Fig. 3(a) (solid lines), reported in 
panel (a) and (b), respectively 
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FIG. 6: (Color online) Semi-log plot of the maximum values, 
within the layered medium bulk, of the normalized squared 
field components as a function of the normalized input field 
amplitude obtained for the full-wave evaluation of the slab 
transmissivity reported in Fig. 5 (a). 



the presence of a TM radiation of free-space wavelength 
A — 810 nm. We have chosen c?i = 2 nm, d2 = 5.25 nm, 
L = 2/ko = 258 nm, ei = -28.81 + lOi, £3 = 10.9 - 3.8i, 
= 3.16 X 10~^^ m^/V^ and X2 = 0. The parameters 
of medium 1 coincide with those of silver |15l |, c harac- 
terized by a very large nonlinear susceptibility [la] , with 
the imaginary part of the permittivity corrected by the 
layer size effect (since di = 2 nm) [T^] whereas medium 2 
is a linear dielectric with gain (to compensate the metal 
losses). The effective permittivity and nonlinear suscep- 
tibility of the considered sample are e = —0.054 -I- 0.007« 
and X — 8.72 x 10~^'^m^ /V^. In Fig. 5 we report the com- 
parison between the transmissivity evaluated through 
full-wave simulations (dotted line) and the transmissivity 
obtained by solving Eqs.([2]) (solid line) in the two situa- 
tions (panel (a) and (b)) corresponding to those of Fig. 2 
and 3. We note that good agreement exists between the 
results of the two kind of simulations and, most impor- 



tantly, that finite-element simulations still predicts the 
above discussed multistability thus proving its robust- 
ness. The origin of the discrepancies between the trans- 
missivities lies in the fact that a surface plasmon reso- 
nance occurs when the TM wave impinges on the slab 
as reported in Fig. 1(b). Therefore an electric field y- 
component (Ey) of plasmonic origin arises mainly at the 
edges of the layers inside the medium and it is character- 
ized by a sub- wavelength varying profile and an evanes- 
cent field tail in vacuum which does not contribute to 
the power flow. Evidently, if Ey is much smaller than E^ 
and Ez, the homogenization theory correctly describes 
the slab nonlinear behavior. In Fig. 6 we report the max- 
imum values (within the medium) attained by the three 
field components for 6 = 0.17 rad as a function of the 
normalized input field amplitude corresponding to the 
evaluated full-wave transmissivities reported in Fig. 5(a). 
It is evident that, outside the shaded region, Ey (dot- 
dashed line) is much smaller than both E^ and E^ (solid 
and dashed lines respectively) and in the corresponding 
regions of Fig. 4(a) the agreement with the homogeniza- 
tion approach is very satisfactory. Within the shaded 
region of Fig. 6, Ey is not negligible and this partially 
breaks the validity of the homogenization approach thus 
leading to the discrepancies of the transmissivities in the 
range 0.02 < |x|^/^|i?j| < 0.04. reported in Fig. 5(a). 

The intensity of the incident plane wave is / = 
{I /2) eo/ fJ.a\Ei\'^ which, for the amplitude range 0.02 < 
x|^/^|i?i| < 0.04 where multistability occurs for e = 
—0.05 (see Figs. 2(a) and 4(a)) and for the above ef- 
fective nonlinear susceptibility X: yields the values 
1.2 MW/cm'^ < / < 4.8 MW/cm'^, which are intensities 
smaller than those normally required for observing the 
standard optical bistability. However it is evident that 
the smaller |e| the smaller the intensity required for ob- 
serving the two-parameter hysteresis so that very smaller 
intensity thresholds (of the order of W/cm^) are very 
likely to be attained. 

In conclusion, we have shown that a nonlinear meta- 
material slab with a very small effective linear permit- 
tivity can host more than one electromagnetic configu- 
ration compatible with the incoming field and this leads 
to multistability and two-parameters hysteresis steered 
through the input intensity and incident angle. The two 
degrees of freedom potentially offers an infinite number 
of hysteresis loops. In addition, as a consequence of the 
small value of the permittivity, hysteresis is predicted at 
low input intensities and this suggests that our findings 
can have important applications in devising tunable and 
compact photonic devices for advanced optical control 
and memory functionalities. 
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